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^— ^ Einstein field equations under spherically symmetric space-times are considered here in 

connection to dark energy investigation. A set of solutions are obtained for a kinemat- 
ical A model, viz., A ~ (a/a) 2 without assuming any a priori value for the curvature 
constant and the equation of state parameter u). Some interesting results, such as the 
' nature of cosmic density f2 and deceleration parameter q, have been obtained with the 

5H consideration of two-fluid structure instead of usual uni-fluid cosmological model. 

Keywords: dark energy, A-model, curvature constant. 

PACS Nos.: 98.80.Jk, 98.80.Cq. 



1. Introduction 

Present-day cosmological research hovers around the investigation of dark energy, an 
exotic type of entity responsible for genera ting acceleratio n in expanding Universe. 
In fact, various recent observational results Hl 2 | 3 | 4 | 5 | 6 H|£| SU gg es ^ t hat the Universe 
is expanding with an acceleration alone while some other works indicate 
that the acceleration is a phenomenon of recent past and was preceded by a phase 
of deceleration. 

Now, the exact nature of dark energy being still unknown, its investigation are 
going on along various paths. Phenomenological models are also contenders in this 
dark energy investigation. Although these type of phenomenological models do not 
originate from any underlying quantum field theory, yet they are useful enough to 
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arrive at some fruitful conclusions. Out of three main variants of phenomenological 
models, viz. kinematic, hydrodynamic and field-theoretic models ^21, the present 
work deals with kinematical models where the dar k en ergy representative A is as- 
sumed to be a function of time. Recently, Ray et alff^ and Mukhopadhyay et alP^I 
have shown the equivalence of four A models, viz. A ~ (a/a) 2 , A ~ a/a, A ~ p 
and A ~ H for spatially flat (k = 0) Universe. But, since the closed (k = 1) and 
open (k = — 1) Universes cannot be entirely ruled out, so there is enough reason to 
investigate dark energy for general k. In this work, therefore, one of the equivalent 
A models, viz. A ~ (a/a) 2 is selected to solve Einstein equations for general k in 
order to have a broader view of accelerating Universe. 

The scheme of the investigation is as follows: Sec. 2 and 3 deals respectively with 
the Field equations and their solutions while some physical feature arising out of 
this work are described in Sec. 4. Finally, some conclusions are made in Sec. 5. 



2. Field Equations for the Spherically Symmetric Space-times 

The Einstein field equations are given by 



R' 



1 



Rg* 



-8nG 



A 
8^G 
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where the cosmological term A is time-dependent, i.e. A = A(<). 
Let us choose the spherically symmetric FLRW metric 

dr 2 



ds z 



-dr + a{ty 



l {d6 z 



(1) 



(2) 



1 — kr 2 

where k is the curvature constant and a = a(t) is the scale factor. For the metric 
given by equation (2), the field equations (1) yield Friedmann and Raychaudhuri 
equations respectively given by 



k 8nGp 
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(3) 



(4) 



where c, the velocity of light in vacuum, is assumed to be unity. 

The generali zed energy conservation law, when both A and G vary, is derived 
by Shapiro et alP^ using Renormalization Group Theory as well as by Vereschagin 
and YegoriarffSl using Gurzadyan-Xue formula^. Since in the present work G is 
assumed as a constant and A is a variable, then the above mentioned generalized 
conservation law reduces to the particular form 



8ttG( P + p)[- 



8ttG 



A 



(5) 



The barotropic equation of state relating pressure and density is given by 

p = ujp (6) 
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where the barotropic index u> can assume the values 0, 1/3, 1 and —1 for pressure- 
less dust, electromagnetic radiation, stiff fluid and vacuum fluid respectively. 
From (1), using equation (6), we get, 



A a 



4ttG(1 + 3c) V3 3 



(7) 



Again, differentiating equation (3) and using equations (5) - (7) we obtain the 
differential equation 



a\ 2 /a\ k _ ( l + w 

a) + 1 + 3c U / + a? ~ I i + 3 W 



(8) 



3. Solutions for the Phenomenological Model A = 3a(a/a) 2 

Using the ansatz A = 3a(d/a) 2 , we immediately get from equation (8) 

- = [3a(l + c) - (3c + 1)] f - (3c + 1) A. 
a la Zaa 

The above equation after simplification reduces to the form 



. d 
aa— 

dt 



ln(da- s / 2 ) 



(3c + l)fc 



(9) 



(10) 



where s = 3a(l + c) — (3c + 1). 

Let us now study the following case when s = —2. In this case equation (10) 
reduces to 



. d (3c + l)fc 
aa dt " n( - aa " = 2 ■ 

Solving equation (11) we get our solution set as 



a(t) - 



C 't + Ci' - 



(11) 



(12) 



H(t) 



Co' - (1 + 3c)fci 



C 't + Ci' - ^±i)fci 2 



(13) 



P(t) = 



3(1 - 3a) 



c ' - i*£ii*t 



2fc 



Co'* + Ci 



/ _ (3^+i; 
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A(t) 



3a[C ' - (3c + l)kt} 2 



4[Co't + Ci' - ^±11 fci 2 ] 2 
where Co' = 2Cq, Ci' = 2Ci, Cq and Ci being constants of integration. 
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If we impose the boundary condition a(t) =0 when t = 0, then C\ = which 
implies C\ = 0. Then the simplified solution set becomes 

1/2 

a(t) 



H(t) 



Co' - (1 + Zu)kt 



C 't - ^^fci 2 



(16) 
(17) 



p(t) 



, [(a + 1)Cq' - (1 + 3w)fci 2 + 2(1 + 3uj)kCo't - ^§^-kt 2 } 
16ttG(1 + 3u)[C 't - ^§^-kt 2 } 2 



(18) 



A(t) = 



3a [C ' - (3w + l)kt]' 



C 't - ^tii fc< 2 



(19) 



It is clear from equation (19) that for a repulsive A, a must be positive whereas 
a = implies a null A. This means that we are getting Einstein's expanding Universe 
without A. 



4. Physical Features of the Solutions 
4.1. Density of the Universe fi 

The above solution set is obtained by assuming s = — 2. Now, s = — 2 means 

2 1 

(20) 



3(l-a)(l+w) 2' 

For k = we get from equations (16), (17), (19) and (20) respectively a(t) oc 
f 2/3(i- Q )(i+u)o-( t ) ^ 1 ^ fj ^) K and oc 1/i 2 . These results were obtained 

by Ray et all^ for flat [k = 0) Universe. Again, from equation (20) we have 

7^ = 3 - (21 » 

Since equation (19) suggests that for a repulsive A we must have a > 0, then 
from equation (21) we find that either u> > 1/3 or u < — 1. For u > 1/3 we get a 
Universe where contribution of electromagnetic radiation is negligible (for radiation 
dominated Universe, w = 1/3) while u> < — 1 signifies the presence of phantom 
energy. Again, using equation (18), the expression for cosmic matter energy density 
Q m can be easily derived and is given by 

n 2{a + 1) i ^-^^ (22) 

Um ~ (1 + 3c) + [Co' - (1 + 3cj)kt} 2 ■ 
Also, from the ansatz A = 3a(d/a) 2 we get 

Q A = a. (23) 
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Then, using equations (21) 



- - r ( 



n A = i 



(23) we obtain 



(24) 



[Co' -(l + 3uj)kt] 2 ' 

For the flat Universe (k = 0), equation (24) reduces to the case of Ray et alP^l. 
Again, equation (24) shows that at time t = 0, the sum of Q m and £Ia becomes 
independent of the curvature constant k and takes a unit value whatever may be 
the value of k. On the other hand, when t tends to infinity, from (24) we have 

n m + n A = i-—?— , (25) 

1 + 6UJ 

From equation (25) we again observe that, f2 m + fiA is independent of k. Thus, both 
the early and late phases of the Universe exhibit the same behaviour so far as the 
curvature dependency of the sum of fl m and f2 A is concerned. It has already been 
shown that for physically valid a, either u) > 1/3 or u> < —1. In the former case, 
2/(1 + 3u>) < 1 and hence by equation (25) 

< fi m + A < 1. (26) 

But, for u> < —1 we have —2/ (1 + 3uj) < 1 and hence equation (25) provides the 
following constraint 



1 < Sl r 



ft A < 2. 



(27) 



The above two relations (26) and (27) suggest that in distant future not only 
the sum total of matter and dark energy density will be independent of curvature 
of space but also they will be either less than (for uj > 1 /3) or greater than (for 
to < — 1) unity which misfits with the present status of the sum of two type of 
energy densities dominating the present Universe . This result is very important for 
visualizing the cosmic evolution in future. 

Now, let us suppose that the Universe is composed of a mixture of two types of 
fluids having barotropic indices uj a and cob (say). Then, from equation (25) we get 

(O m + fl A ) a = 1 - 2 , (28) 
1 + 6uj a 



(fl m + Q\)b = 1 



1 + 3u>b 



If (fl m + fl\)avg be the average of (fl m + fl A )a and (fl m + fl\)b then 



(i7 m -\- Q\*)avg 1 



2 + 3(uj a + Lo b ) 



Xl + 'iuo a ){l + 'iuob) 
Equation (30) shows that when uj a + LOb = —2/3, then 

(^m ~t~ ^A^)avg 1- 



(29) 



(30) 



(31) 



This means that like the early and present Universe, for late Universe also the sum 
of fi m and J7a will be unity only if the Universe contains a mixture of two types 
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of fluids rather than a single fluid. Since uj a + ujb = —2/3 and it has already been 
shown that either uo > 1/3 or w < —1, then let us suppose, u a = 1/3 + e and 
lui, = — 1 — e where e > 0. Therefore, for large value of t, when the Universe is filled 
with a mixture of two types of fluids (one of them being phantom-fluid) and if the 
value of barotropic index of one fluid is — (1/3 + e) and that of the other is (— 1 — e) , 
then the average value of the sum of fl m and Q\ will be equal to one. 

4.2. Deceleration parameter q 

Let us now consider the expression for the deceleration parameter q which is given 

by 

«~S~(> + £)- < 32 > 

If the Universe is composed of two fluids with equation of state parameters uj a 
and cjf,, then each of them will have some effect on the dynamics of the Universe. 
So for calculating the value of the deceleration parameter q, contributions coming 
from each component should be taken into account. If q a and be the values of 
two separate parts of q coming from fluids having barotropic indices uj a and uii, 
respectively, then using equation (17) we get from equation (33) the following two 
expressions 

Co' 2 

Qa = [Co'-(2 + 3e)fct] 2 ' (33) 

C ' 2 

Qb = [Co' + (2°+3e)fcip- (34) 

If q e ff be the effective value for q, coming after considering the separate parts q a 
and qb, then 

4Co ,3 fcf 

Qeff = To • ( 35 ) 

[Co' - (2 + 3e) 2 fc 2 i 2 ] 2 

Equation (35) shows that the sign of q depends only on two quantities, viz., the 
integration constant C/(= 2Co) and the curvature constant k. If for simplicity we 
assume Co to be positive then we get an accelerating or a decelerating Universe 
according as k < or k > 0. This result can be interpreted as follows. The Universe 
is made of two types of fluids having equation of state parameters u a and w&, one of 
which is acting as a prohibitor and another as a supporter of cosmic acceleration. 
In the previous matter dominated phase, k had a small positive value (i.e. q a > q^) 
and as a result the Universe was decelerating. But at a certain time during cosmic 
evolution, the second type of fluid (viz., phantom fluid) took the upper-hand (i.e. 
Qa < Qb) and consequently k has become slightly negative. That is why the present 
Universe is in a state of acceleration. A very small positive or negative value of 
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the curvature constant do not contradict the observational result that the present 
Universe is nearly flat. Further, the change of sign of q shows that the cosmic 
acceleration is a recent phenomenon. 

5. Conclusions 

The present work, apart from being a generalization of an earlier one^l, has revealed 
some new and interesting physical features also. For a flat Universe, all the results of 
Ray et alff^can be recovered from the expressions of a(t), Hit), pit) and A(t) of the 
present work. Moreover, it has been possible to trace the entire cosmic evolution, 
starting from the Big-Bang and extending to distant future. The most significant 
result is related to the cosmic matter and dark energy densities for non-flat Universe. 
It has been shown that for non-flat Universe, f2 m + Qa = 1 only when the Universe 
is composed of two types of fluids, one with u> > 1/3 and another with u> < — 1. 
This means that both stiff-fluid (lo — 1) and phantom- fluid {lj = —1) can be one of 
the two constituents of the Universe when the Universe is not flat. The evolution of 
the deceleration parameter q shows that a non-flat Universe would be decelerating 
in the past and accelerating at present. This result is very significant for A-CDM 
cosmology. 

Another interesting poi nt is the absence of Big-reap even in the presence of a 
fluid with it) < — 1. Caldwelff^ and Caldwell et alP2 demonstrated the occurrence of 
a Big-reap in the presence of a fluid with supernegative equation of state. It may be 
mentioned here that Gonzalez-Diaz h as shown that by a proper generalization 
of the Chaplygin-gas model, a Big-r eap can be avoided even in the presence of 
phantom energy whereas Abdalla et al)^' have arrived at the same result through a 
slight modification of GTR. But in the present work, a cosmic doomsday is shown 
to be impossible within the normal framework of GTR. One of the reasons of it may 
be the presence of another fluid apart from the phantom fluid. Th at o ther fluid may 
act as an inhibitor of Big-reap. It may be mentioned that Stefancic^^has developed 
a model in which the dark energy component and the matter component interact 
with each other resulting in the appearance of phantom energy out of non-phantom 
matter. The present work can be considered as a counter example of that because 
here phantom matter, in the presence of another component as a mediator, can 
behave as a non-phantom matter. 

Finally, it is to be noted that the present work has demonstrated that although 
current observational data points towards a k = Universe, yet we are not in a 
position right now to completely rule out k = ±1 cosmologies. 
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